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ABSTRACT 

Context. In support of NASA's exploration program and the return to the Moon, the polarimetric signature of dispersed individual 
Lunar regolith dust grains is studied to enable the characterization of the dust exopsheric environment by remote, in-situ, and standoff 
sensing. 

Aims. We explore the value of the negative polarization branch (NPB) as a signature for characterizing individual grains to determine 
if it can be used in the same way as for surfaces of planets and atmosphereless bodies. 

Methods. The linear polarization phase curve for single spheres of silicate and for aggregates of spherical silicate grains of different 
porosity are computed using the discrete dipole approximation (DDA) for a range of grain sizes. Features in these curves are identified 
and their evolution explored as a function of grain size and aggregate porosity. We focus particularly on the so-called negative 
polarization branch that has been historically used to characterize planetary surfaces. 

Results. Calculations show that polarization phase curves for spherical grains exhibit a sharp transition over a narrow range of 
size parameter between two distinct regimes, one typical of Rayleigh scattering and another dominated by a large NPB. The linear 
polarimetric signature observed for aggregates is a composite of a) the polarization induced by individual grains composing the 
aggregate and b) the polarization due to the aggregate as a whole dust grain. The weight of each component varies depending on 
the porosity of the aggregate. An NPB similar to the one observed for atmosphereless astronomical bodies is present for different 
ranges of the size parameter depending on the value of the porosity. It appears as a remnant of the negative branch exhibited by the 
single spherical grains. The sharper, narrow negative branch that is measured for some granular surfaces in the laboratory or seen in 
astronomical observations is not observed here. 

Conclusions. These results suggest that the wide negative branch is due to the scattering by individual grains and single aggregates, 
while the narrow negative branch is more likely due to coherent backscattering or shadowing effects in bulk material. The shape and 
evolution of the NPB could be used to characterize spherical grains and to differentiate between aggregates with the same porosity but 
different sizes, but does not appear to be a practical candidate for univocally differentiating between aggregates of different porosity. 
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1. Introduction 

The nature of the Lunar dust environment is an unknown that 
needs to be understood before long-term human exploration and 
settlement can become a reality. A return to the Moon will in- 
evitably be accompanied by anthropogenic transport and dy- 
namics. The transcripts and debriefings of the Apollo missions 
demonstrate the pervasive nature of lunar dust in an exploration 
environment: invasion of the Lunar Module, solar panel and ra- 
diators clogging, and the possible induction of allergic reactions. 
Eugene Cernan, Mission Commander of Apollo 17 summarized 
the issue : "The dust was very difficult to work in and was a 
big hindrance. It obscured your vision if you let it get on your 
visor. It was almost like a dew or mist because it clung to ev- 
erything. It got on every movable surface, it got in the suits, and 
when we got out of our suits in the spacecraft in between EVAs, 
it got in the pores of our skin and got under our fingernails. 
And it didn 't just get on the outside parts of our nails and get 
them dirty but, literally, it got down between the skin and the 
nail. It took three months for lunar dust to grow out from under 
my nails. It infiltrates." (Apollo Lunar Surface Journal dALSJI ), 



Apollo 17 Journal). The composition and morphology of lunar 
dust have made it suspect of being toxic to humans, and its prob- 
able transport into future habitat modules and inhalation by the 
crew could prove to be more than a simple inconvenience. Other 
physical characteristics (such as high abrasiveness) are respon- 
sible for its position near the top of the list of engineering chal- 
lenges. NASA has acknowledged the importance of these issues 
by forming the Lunar Airborne Dust Toxicity Advisory Group 
(LADTAG), which regroups interdisciplinary experts and has 
been charged with setting health standards and risk criteria for 
use by engineers and astronauts during lunar mission design and 
operation. 

In addition to the transport generated by human activity, it 
has been suggested that dust grains in the micron and sub-micron 
size range can be naturally transported by various electrostatic 
mechanisms above the surface of the Moon to altitudes in excess 
of 100km (e.g., Stubbs et al. 12007b . Evidence of a dust exop- 
shere includes observations of "streamers" and "horizon glow" 
by Apollo Astronauts (McCoy & Criswell 1974) as well as by 
the Surveyor landers (Rennilson and Criswell [l9741 l and Apollo 
coronal photography (McCov fl976l ). The scale-height of the ob- 
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served horizon glow has been estimated to be ~10km (Zook & 
McCoy [1992 Murphy & Vondrak [19931 The dusty lunar sur- 
face is electrostatically charged by the photoemission of elec- 
trons by solar ultraviolet, as well as the incident electrons and 
ions from the ambient plasma environment. The resulting sur- 
face electric fields drive the transport of charged lunar dust. 
Different mechanisms have been proposed, from the levitation 
of micron-size grains at ~30cm above the surface, to the lofting 
of submicron grains to high altitudes by the dynamic fountain 
effect (Stubbs et al. 2006). These fundamental transport mech- 
anisms are also expected to be active on other atmosphereless 
solar system bodies with granular surfaces, making the results 
presented here relevant to other objects such as asteroids and the 
surfaces of airless planets and planetary satellites. 

The capability to characterize the dust environment near the 
surface and in the atmosphere of the human habitat is critical to 
enable a human settlement; characterization of the grain popula- 
tion lofted to higher altitudes is equally important for our under- 
standing of the fundamental transport mechanisms. Our current 
work aims at developing technology to probe and characterize 
this environment at all altitudes by in-situ, standoff and remote 
sensing; the results reported in this publication are a preliminary 
step towards building a knowledge base to support that objective. 

While early polarimetric observations of the Moon were at- 
tempted by Arago ( 1181 II |1844| l, the doctorate thesis of Bernard 
Lyot d 19291 1 laid the foundations for comparative studies of the 
polarimetric properties of celestial bodies and terrestrial sub- 
stances. Measuring the optical properties of various materials- 
including ash from several eruptions of mount Vesuvius-and ob- 
serving the light scattered by astronomical objects he noticed 
that the polarization curves of the Moon, Mercury, certain ar- 
eas on Mars and asteroids most resemble the curves obtained 
for opaque powdered materials in the laboratory. Supported by 
observations of integrated moonlight and regional studies, he 
concluded that the Moon is covered almost entirely by a pow- 
dery material closely resembling terrestrial volcanic ash. One 
particular feature leading to this conclusion was the presence of 
a negative branch on their polarization curve at small phase an- 
gles (The sign of P depending on which linear polarization com- 
ponent dominates the observed reflectance.) His work was fol- 
lowed by that of Dollfus ( fT955irT998irT999l[2000l >. who notably 
selected by polarimetry a fine grained, very cohesive basaltic 
material that reproduced satisfactorily the astronomical obser- 
vations (Dollfus [19661 ); polarization curves from the several ar- 
eas and features of the lunar surface were described in terms 
of specific parameters (such as the minimum and maximum of 
polarization, the inversion angle, the gradient at the inversion an- 
gle) and compared to laboratory measurements on meteoritic and 
mineralogic samples. The first samples from the lunar surface 
returned by Apollo 1 1 showed a striking resemblance with the 
selected material (Dollfus [19981 ), thus validating the technique. 
{Numerous other works have generally focused on inferring the 
properties of particulate surfaces through the characteristics of 
their observed polarization phase-curves; Geake & Geake ( 1990) 
and Kolokolova et al. (1993 ) studied how the polarization-phase 
curve of powder surfaces varies in scale and shape with grain 
size for sub-wavelength grains; Shkuratov et al. (2002 2007a, 
2007b) measured the scattering by particulate surfaces of vari- 
ous compositions and focused on the origin of the opposition ef- 
fect and negative branches of the phase curve; Lasue et al. ( 2006 ) 
and Hadamcik et al. (2003 , 2003 ) modeled cometary light obser- 
vations by the light scattered by porous irregular aggregates of 
sub-micron grain sizes of silica and organics. 



The origin of the negative branch, commonly associated with 
powdered materials and fluffy aggregates, is still debated. Lyot 
proposed that negative polarization is due to three mechanisms : 
multiple reflection, refraction in transparent particles and diffrac- 
tion. While these are very broad and general, more precisely de- 
scribed mechanisms include shadowing or coherent backscatter- 
ing in fluffy grains. The Mie solution to Maxwell's equations- 
Mie theory-for samples of polydispersed separated spherical 
grains of different refractive indices have been proposed to ex- 
plain the negative polarization of cometary atmospheres. These 
models have not successfully been applied to planetary regolith. 
Shkuratov et al. ( 1 19941 ) provides a detailed review of these var- 
ious models and mechanisms. In addition to the wide negative 
branch observed by Lyot, a narrow peak (less or of the order of 
a degree angle wide) of negative polarization is sometimes ob- 
served in conjunction with a sharp increase in scattered intensity 
(opposition effect) at angles close to backscattering. It is still un- 
clear whether both the wide and narrow negative branches have a 
common origin or are due to different mechanisms. Ovcharenko 
et al. ( 120061 ) proposes-based on laboratory measurements of 
granular surfaces-that the wide negative polarization branch is 
due to individual aggregate scattering while the narrow branch 
is due to coherent backscattering on surface layers. 

Scanning electron photomicrographs of Apollo samples re- 
vealed the presence of a variety of grain morphology, from mi- 
crometer and sub-micrometer agglutinates with irregular and 
sharp edges-which origin lies in the surface bombardment by 
meteorites-to smoother glass droplets of volcanic origins (see 
McKay et al. (1991} for a review of the lunar regolith proper- 
ties.) But it wasn't until recently that the fine (from lOOnm to 
20jt/m) and ultrafine (20 to lOOnm) particle content of the lu- 
nar regolith was measured (Greenberg et al. 2007 ), the difficulty 
residing in the separation and counting technique that need to 
be used for sub-micron particles. Along with measuring Particle 
Size Distribution functions (PSD), Greenberg and collaborators 
also took hundreds of photomicrographs of grains; typical par- 
ticles in the submicron regime-the most relevant for electrostat- 
ically transported grains-do not display any complex feature or 
any departure between aerodynamic and physical dimensions. 
They are compact ellipsoid particles with small aspect ratios. 
In only two photomicrographs did they observe aggregates of 
these compact particles. It is unknown though what proportion 
exists in aggregate form in the lunar environment, as the protocol 
to obtain PSD functions included a de-agglomeration process. 
Spherical grains and their aggregates are therefore a relevant ap- 
proximation to electrostatically transported dust. 

The present work initiates a systematic study of the polari- 
metric properties of individual grain models in order to charac- 
terize the lunar exospheric environment, where individual dust 
grains are dispersed above the surface. It is therefore different 
in nature than bulk or surface regolith characterization. We first 
describe the computation method and target models; we then 
present the linear polarization due to the scattering by single 
spheres and follow with the results for random spherical grain 
aggregates of various porosity. For each target type, a range of 
sizes is explored. We discuss the origin of the NPB and its use- 
fulness as a dust grain characterization tool. 

2. Model and Approach 

Computations have been carried with the publicly available code 
DDSCAT developed by Draine & Flatau (2004) implementing 
the Discrete Dipole Approximation (DDA) first proposed by 
Purcell & Pennypacker (119731 ). The DDA approximates a con- 
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Fig. 1. Sensitivity of the model to variations of the complex re- 
fractive index. 

tinuum target by a finite array of polarizable points, therefore 
allowing to calculate the scattering and absorption properties of 
targets of arbitrary geometry. Each point acquires a dipole mo- 
ment in response to the local electric field induced by the in- 
cident light wave. These dipoles interact with one another via 
their electric field-which is why the DDA is also referred to as 
the coupled dipole approximation. The method is described in 
more details and reviewed in Draine & Flatau ( 119941 1. The major 
inputs to the code are the complex refractive index, the polariza- 
tion state and frequency of the incident wave, the size, geometry 
and orientation of the target; outputs include the 4x4 Mueller 
scattering intensity matrix S , ; , as well as scattering, absorption 
and extinction efficiencies. The linear polarization degree is de- 
fined in terms of the Mueller matrix elements by P = -Sn/Su- 

While the signatures of individual grains of the various 
shapes and compositions known for lunar material will eventu- 
ally have to be determined in the future, we restrict the subject of 
this publication to single spherical grains and aggregates of 28 
spherical grains. For aggregates, we define the porosity param- 
eter as the percentage of void space in the target defined by the 
aggregate inscribed in the smallest possible box. Computations 
were carried out for three values of this parameter: 0.9,0.8, and 
0.6. To the authors knowledge, the only measures of refrac- 
tive index of lunar regolith are these of plagioclase grains from 
Apollo samples, of size ranging from 0.75mm to 1 .55mm. These 
grains have a real refractive index between approximately 1 .56 
and 1.59 (King et al. 1 1 97 li t. Because the size of these samples 
is much larger than the fine grains that interest us, we chose a 
complex refractive index that reflects the optical properties of 
typical astronomical silicate dust, m = 1.62 + i0.003 (Lasue & 
Levasseur-Regourd 2006). The sensitivity to both the real and 
imaginary part of the refractive index can be seen on Fig. Q] 
where the polarization phase curve for an aggregate of poros- 
ity 0.6 is shown for various values of m. This demonstrates that 
practical applications of polarimetric signature characterization 
of the lunar dust environment will require a better knowledge 
than we have today of the optical properties of the fine and ultra- 
fine portion of the regolith. With the recent renewed interest in 
Lunar Science, we do not expect for this gap in our knowledge 
to persist for long. 

The size parameter is varied for each kind of target; it is de- 
fined as 
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Fig. 2. Notations used in this paper to describe polarization curve 
characteristics. 



where k is the wavenumber and A is the wavelength of the in- 
cident radiation, and a is the radius of the sphere of equivalent 
volume. For aggregates, we define X agg as the size parameter 
of the whole aggregate, and X, as the size parameter of its con- 
stitutive individual spherical grains. The incident light is unpo- 
larized and a reference wavelength /lo=628.3nm is chosen for 
practicality (X = 2na/Ao = 10a with a in microns). For non- 
spherical targets, scattering is computed and averaged over more 
than a thousand target orientations and is therefore representa- 
tive of single scattering by a sample of separated randomly ori- 
ented grains. For each computation, the validity conditions for 
the DDA model are verified : \m\kd < 1 and d is smaller than any 
structural length scale of the target, where d is the inter-dipole 
separation. We limit the exploration in this paper to size parame- 
ters for which a single negative branch exists in the polarization 
phase curve. We use the scattering angle rather than the phase an- 
gle (which is more customary in astronomy.) These two angles 
are supplementary; a scattering angle of 0° corresponds to for- 
ward scattering and 180° to backscattering. P„„„ and P mflA refer 
to the maximum and minimum value of the polarization curve, 
and m i„ and 9 max the scattering angles for which these values 
are reached. The inversion angle 0,- m , is the angle for which the 
polarization curve changes sign. These notations are illustrated 
on Fig. [2] 

The present work is concerned with spheres and sphere ag- 
gregates of homogeneous composition; the refractive index is in- 
dependent of wavelength. More complex shapes, compositions, 
spectral properties as well as different aspects of the polarimetric 
signature (beyond the characteristics of the negative polarization 
branch) will be considered in future publications. 

3. Polarization by single spheres 

The linear polarization phase curves for spherical grains of four 
representative size parameter are shown on Fig. [3] For small 
values of X (for example X = 1.0 in this figure, correspond- 
ing at this wavelength to a grain of size 0.1 fmi) the phase curve 
is expectedly similar to that for Rayleigh scattering (for scat- 
terers of size smaller than the wavelength): a bell-shaped curve 
peaking close to 90° scattering angle. As the size parameter in- 
creases, higher modes develop; a negative branch appears around 
X — 1.7; forX = 4.0, the curve exhibits four local minima, and 
for X = 10.0 it shows 8. This behavior is classically known and 
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Fig. 3. Linear polarization (P) of light scattered by silicate 
spheres of variable size parameter (X). 

is theoretically described by the Mie theory (Hulst [l98Tl ). As the 
behavior of the polarization curve becomes more complex, sev- 
eral 'negative branches' appear. In the present article, we will 
focus exclusively on scattering regimes where no more than one 
negative polarization branch (NPB) is present. For the present 
computations, this translates to a size parameter X < 2.5. 

Figure |4]explores the polarimetric signature and particularly 
the characteristics of the NPB in this size range. Figure Ufa) 
shows the evolution of the shape of the phase curve through five 
representative values of X as it continuously transforms from 
Rayleigh-typical to a curve dominated by a negative branch that 
is both wide and deep (with Q mv close to 130° and P m ,„ close to 
-1 at X — 2.) The Rayleigh-typical regime extends roughly until 
X = 1.5; the large NPB regime is observed for X > 1.8. This is 
clearly illustrated by Fig.Hfb) showing the minimum and maxi- 
mum of the polarization phase curve P„„„ and P ma x as a function 
of X; for X < 1.5, P m ,„ is equal to zero and P max is close to 1; 
for X > 1.8 P max averages a much lower value of 0.15 and P„„„ 
reaches close to - 1 . 

As a precursor to the transition between the two regimes, 
the scattering angle for which the polarization reaches its max- 
imum, 9 max , first increases up to 120° (FigHJc)). The switch is 
then abrupt, reminiscent of a phase transition: the polarization 
maximum P max starts to collapse for grain sizes close to X = 1 .57 
(marked by a vertical dotted line on Figj4fb) and (c)); 8 max then 
reaches a maximum (at X — 1 .62) and the curve peak starts to 
shift to smaller angles (Figj4{c)); when P max has reached half its 
initial value, the negative branch appears (forX - 1.65, marked 
by the dashed lined of FigUb), (c) and (d)) and rapidly grows 
such that for X = 1.8 it has reached down to P„„„=-0.95 at an 
angle Q mm of 130° (FigHfc)), with an inversion angle 6>, m , of 95° 
( FigHd)). 

For values of X > 2.2 the appearance of additional features 
of the curve raises the value of P„„„ (see the curve for X = 2.4 
on FigHfa) for illustration). Polarization by larger grains shows 
very different signatures (as shown on PigOJ an d will be the ob- 
ject of future work. 

This transition over a narrow range of size parameter X = 
2na/A could provide a valuable signature for the characteriza- 
tion of this type of grain for size and composition. For example, 
by rotating scanning frequencies and finding the frequency for 
which this transition occurs, composition and/or size of spheri- 
cal grains could be deduced. 
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Fig. 4. Linear polarization (P) phase curve of light scattered by 
spheres of silicate of variable size parameter (X); (a) P as a func- 
tion of scattering angle, for five values of X; (b) Polarization 
maximum P ma x and minimum P m ,„ as a function of X; (c) Angles 
for which the maximum (f?„ ra v ) and minimum (#„„■„) are reached 
as a function of X; (d) Inversion angle 8j m ,, indicating the width 
of the Negative Polarization Branch as a function of X. Vertical 
dashed lines mark the value of X for which the NPB appears. 
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Fig. 5. Prototype aggregates of porosity 0.9 (left) and 0.8 (mid- 
dle) and 0.6 (right). Constitutive grains are of equal size. 

4. Polarization by aggregates 

Lunar regolith is partly composed of aggregates of various size, 
density and composition. As a first step towards a realistic de- 
scription of these dust grains, we have computed the polariza- 
tion phase curves for random aggregates of spherical grains. 
Aggregate prototypes have been created of three different poros- 
ity : 0.9, 0.8 and 0.6. These aggregates are shown in Fig. [5] As 
in the single sphere case, scattering is computed for each target 
type for a range of size parameters. Individual grain size scales 
with the aggregate. (For a given porosity, if aggregate B is twice 
smaller than aggregate A, the grains in B are also twice smaller 
than in A.) This choice is voluntary and does not reflect a limita- 
tion of the model or the DDA code. In this light, various values 
of X for a given aggregate porosity can be seen as experiments 
where 

- either the incident wavelength is constant and illuminates 
different aggregates composed of individual grains of vari- 
able size; 

- or the physical size of the aggregate is constant and the inci- 
dent wavelength is variable. 

The linear polarization phase curves for the three porosity 
values are shown in Fig. [6] for several values of the constitu- 
tive grain size parameter Xj. As the size parameter increases as 
well as when the porosity decreases, the general trend is (1) the 
appearance of higher modes, leading to an increased irregular- 
ity and (2) a decrease in the amplitude of the polarization de- 
gree curve. For Xj = 1 and for a porosity of 0.9 the curve very 
much resembles Rayleigh scattering; As the porosity decreases, 
the maximum shifts to larger scattering angles, a secondary peak 
appears between 40° and 50° and the curve becomes bi-modal. 
Further evolution can be seen on the curve for aggregate porosity 
of 0.6, at X/ = 1.29, where a third mode is apparent (Fig. [6jc)). 

The regime dominated by a large and deep NPB observed 
for spherical grains is replaced for aggregates by shallower, ir- 
regular polarization curves; the width of the NPBs remains com- 
parable to the spherical case. The evolution of the polarization 
phase curve resembles the evolution experienced over a much 
larger range of size parameters by the spherical targets (Figj3]l 
but the aggregate curves develop higher modes faster. The am- 
plitude of polarization is lower than in the spherical case. This is 
probably a result of multiple scattering between individual con- 
stitutive grains. Depolarization of an initially polarized incident 
wave due to multiple scattering is a known phenomenon (see for 
example Mishchenko et al. (2006)); as an unpolarized wave is 
polarized by single scattering on the first constitutive grain it en- 
counters, subsequent scattering occurring before its emergence 
from the aggregate will result in some depolarization. 

This behavior suggests that the polarization induced by these 
aggregates can be described as the result of the superposition of 
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Fig. 6. Linear polarization phase curve for aggregates of silicate 
spheres of porosity (a) p=0.9, (b) p=0.8 and (c) p=0.6. For each 
porosity five values of Xi are shown. 



the scattering due to (1) each individual constitutive grain and 
(2) the scattering due to the larger sized aggregate. This is illus- 
trated by Fig|7]showing the polarization curve for aggregates of 
size parameter X agg = 3.9 whose individual constitutive spher- 
ical grains have a size parameter Xj = 1.3 (dashed lines). The 
curves for spheres of size X = 3.9 and X = 1.3 are also shown 
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Fig. 7. Linear polarization phase curves for aggregates of three 
different porosity, of global size parameter X a gg = 3.9 and in- 
dividual grain size parameter Xj = 1.3 (dashed lines), and for 
spheres of size parameters X = 1.3 and X = 3.9 (dotted lines). 
Solid lines are best fit assuming the superposition of single scat- 
tering by a target of size parameter X — 3.9 and multiple- 
scattering up to the ninth order by targets of X = 1.3. 



(dotted lines). For each target, the solid line curve is the best fit 
of the superposition of (1) single scattering by a spherical tar- 
get of X — 3.9, (2) single scattering by spherical targets of size 
X — 1.3, (3) multiple-scattering up to the ninth order by spheri- 
cal targets of size X = 1.3, as described by: 



<P(6) = A P(6,X = 3.9) + Y J A i P n (6,X = 1.3), 



where V is the fitting function, P(6,X) is the polarization 
curve for a silicate sphere of size parameter X and Ao, ...,A„ are 
fitting parameters; n=9 for this figure. Despite the simplicity of 
this formulation, the resulting curves exhibit features similar to 
the ones seen for the actual aggregate. It is evident from these 
curves that the more compact is the aggregate, the more does the 
polarization curve reflect the features of scattering by a spherical 
target of its size parameter. Conversely and expectedly, the more 
porous the aggregate, the closer is its curve to scattering and 
multiple scattering by spheres the size of each individual grain. 

As with the spherical case, the behavior described above is 
quantified by the parameters shown in Fig. [8] The evolution of 
Pmm with X, illustrates the shallowness of the negative polariza- 
tion branches compared to the spherical case; Fig. Oa) shows 
that the evolution of P max with X is more gradual as the porosity 
decreases. Curves representing #,■„,, in Fig.[8jc) show that there is 
no regular monotonic behavior identifiable for the appearance of 
the negative branch as a function of the porosity ; when p varies 
from 0.9 to 0.8 and 0.6, the size parameter for which the negative 
branch appears goes from 1.19 to 1.28 and 1.0. The same con- 
clusion holds for 9„,i„ and 6 max shown in Fig.[8jb). It is neverthe- 
less still clear from these figures that the higher the porosity, the 
closer the curves are to the case of spherical individual spheres. 



5. Conclusions 

We have compared the degree of linear polarization of light scat- 
tered by spheres of silicate and aggregates having three different 




Fig. 8. Characteristics of the linear polarization phase curves of 
aggregates: (b) P min and P max , (c) min and max , (d) 9 inv , as a 
function of X. In (b) and (c) curves are plotted only for the range 
of X for which a single NPB exist. 



values of porosity as a function of grain size. We have focused 
on the range of size parameter for which a negative branch exists 
at high scattering angle (low phase angle). For spherical grains, 
the regime transition over a narrow range of size parameter at the 
onset of the appearance of the NPB makes it a good candidate 
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for single grain characterization; conducting multi-spectral ob- 
servations of the lunar exosphere and identifying the wavelength 
for which a sharp change in polarization occurs would yield in- 
formation about size and composition of grains. For aggregates 
it has been shown that, for high porosity, the behavior is close to 
that of individual spheres. As the porosity decreases, the curve 
becomes closer to that of a larger sphere of size parameter simi- 
lar to that of the whole aggregate. While the NPB can provide a 
useful signature to differentiate between aggregates of the same 
porosity but of different sizes, it would be difficult to differentiate 
univocally between aggregates of different porosity, thus requir- 
ing complementary measurements (e.g. phase function shape) to 
further constrain the geometry. 

Our results suggest that the wide NPB in these aggregates 
is a remnant of that observed in individual constitutive grains. 
Ovcharenko et al. (2006 ) proposes that the wide NPB of plane- 
tary surfaces finds its origin in the individual aggregates, while 
the narrow NPB-which is unseen in our calculations of single 
grains-is the result of coherent backscattering between aggre- 
gates on a granular surface. Combining our results with these 
from Shkuratov et al. ([2007a) and Ovcharenko et al. (2006]), we 
observe that the wide NPB, which is largest for individual spher- 
ical grains, decreases in amplitude for aggregate and decreases 
further for granular surface, essentially originates from the scat- 
tering (single and multiple) by individual grains constituting the 
aggregates. As the polarization for individual aggregates evolves 
from a Rayleigh typical bell-curve to a marked NPB to finally a 
rather flat, low polarization degree curve, it seems probable that 
the negative branch that is always observed for surfaces such as 
the Moon's is caused by the averaging over a distribution of these 
phase curves. Aggregates of small size parameter contribute to 
the positive maximum; large size parameters average to negligi- 
ble contributions; intermediary size parameters contribute to the 
negative branch. 

We have shown the possibilities and limitations of using the 
NPB as a polarimetric signature for the lunar dust environment 
in the simple case of spherical silicate grains and aggregates of 
spheres. Further calculations and analysis are currently being 
pursued to develop a thorough model of scattering properties of 
dispersed lunar regolith taking into account the various shapes, 
structures and agglomeration states found in its granular matter, 
as well as the range of sizes and the varying elemental composi- 
tion. 
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